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ABSTRACT 

In this paper we compare the properties of three sub-systems of Galactic globular clusters, 
which are defined according to metallicity and horizontal branch morphology. We specifically 
focus on cluster luminosities, structures, surface brightnesses and ellipticities. It is shown that 
the so-called "young" halo (YH) clusters, which are thought to have formed in external satel- 
lite galaxies, exhibit characteristics which are clearly distinct from those of the "old" halo 
(OH) and bulge/disk (BD) clusters, the majority of which are believed to be Galactic natives. 
The properties of the YH objects are, in many respects, similar to those of clusters belonging 
to a number of present day satellite dwarf galaxies. The OH and BD populations have appar- 
ently been strongly modified by destructive tidal forces and shocks in the inner Galaxy. By 
comparing the properties of the three cluster sub-systems, we estimate that the present popu- 
lation of native Galactic clusters may only represent approximately two-thirds of the original 
population. Several clusters with low surface brightnesses are observed to be highly flattened. 
We briefly speculate on the possibility that this ellipticity reflects the intrinsic flattening of 
dark matter mini-halos in which these optically dim clusters might be embedded. Finally, we 
examine the distribution of clusters on the size (logi?h) vs. luminosity (My) plane. Three 
objects are seen to fall well above the sharp upper envelope of the main distribution of clus- 
ters on the size-luminosity plane: oj Centauri, M54, and NGC 2419. All three of these objects 
have previously, and independently, been suggested to be the stripped cores of former dwarf 
galaxies. This suspicion is strengthened by the additional observation that the massive cluster 
Gl in M31, plus a number of the most luminous clusters in NGC 5128 also fall in the same 
region of the log i?h vs. My plane. All of the latter objects have previously been suggested as 
the stripped cores of now defunct dwarf galaxies. 
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1 INTRODUCTION 

Our Galaxy contains ~ 150 globular clusters, located at distances 
between 500 pc and 120 kpc from the Galactic centre. Because they 
are generally bright, and probe such a large range of Galactocen- 
tric radii, these objects constitute important tools for both tracing 
the properties of our Galaxy in the present epoch and for piecing 
together its formation history. They are also interesting objects in 
their own right, since they have the potential to provide information 
on star and cluster formation processes. Furthermore they provide 
useful laboratories for the study of self-gravitating stellar systems. 

As a globular cluster evolves dynamically its core contracts 
and its envelope expands. However, A^-body calculations (e.g., 
Spitzer & Thuan (1972); Lightman & Shapiro (1978); Aarseth & 
Heggie (1998)) show that ijj, - the radius that contains half of the 
cluster stars in projection, generally changes very little over periods 
as long as ten relaxation times. This suggests that R^, like [Fe/H], 
may provide valuable information on physical conditions during 
the era of globular cluster formation. However, an important caveat 



is that metallicity dependent differences of up to 20 per cent can oc- 
cur between the half-light radii and the half-mass radii of clusters 
of differing metallicity (Jordan 2004). The reason for this metallic- 
ity dependence is that the masses of the most luminous surviving 
stars become larger as the cluster metallicity increases. As a re- 
sult the density profile of a metal-rich cluster will be slightly more 
centrally concentrated than that of an otherwise similar metal-poor 
cluster. Nevertheless, the half-light radii of clusters are of consid- 
erable interest because they are stable over many relaxation times. 

The work of Mackey & Gilmore (2004) has recently re- 
emphasized the fact that Galactic globular clusters comprise a very 
inhomogeneous class of objects. It is therefore useful to try to iden- 
tify sub-groups made up of clusters having similar characteristics. 
One of the most successful classification schemes for the Galactic 
globular clusters is that introduced over a decade ago (see e.g., Zinn 
(1993); van den Bergh (1993), and references therein), in which 
clusters are grouped according to both their metallicity and their 
horizontal-branch (HB) morphology. These parameters are closely 
related to a number of other cluster properties such as, for example, 
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Galactocentric radius and age. Such a division of globular clusters 
can therefore provide us with a great deal of information about the 
evolutionary history of the Galaxy. Recently, Mackey & Gilmore 
(2004) updated the original classification scheme to cover the full 
Galactic globular cluster sample for which information is presently 
available. By comparing the different Galactic sub-groups to the 
cluster systems in nearby dwarf galaxies (i.e., the LMC, SMC, 
and Fornax and Sagittarius dwarf spheroidal galaxies) these authors 
confirmed that ~ 30 per cent of the Galactic globular cluster sys- 
tem possesses properties consistent with an origin external to the 
Milky Way. In other words, such clusters are most likely to have 
formed within much smaller galaxies (perhaps similar to the Lo- 
cal Group dwarf spheroidals observed at the present epoch) which 
later merged with the Milky Way system. Furthermore, Mackey & 
Gilmore (2004) showed that classifying Galactic globular clusters 
by metallicity and HB morphology leads quite naturally to a segre- 
gation of clusters by size. These authors used cluster core radius R c 
as their size diagnostic; however, unlike the half-light radius Rh, 
this parameter is quite sensitive to the dynamically evolving state 
of a globular cluster. It is therefore of interest to extend the analy- 
sis of Mackey & Gilmore (2004) to additional parameters, such as 
Rh, which "remember" more about the initial sizes of clusters, and 
hence about physical conditions in the Galaxy and/or in its progen- 
itors, at early epochs. 

In the present paper we shall examine the different Galactic 
globular cluster sub-systems in terms of cluster structures, lumi- 
nosities, and surface brightnesses. Section 2 lists the relevant data 
for each of the Galactic globular clusters and presents the classifi- 
cation of each object according to Mackey & Gilmore (2004). In 
Section 3 we use these data to investigate the collective properties 
of each cluster sub-system, while in Section 4 we discuss the possi- 
ble links between half-light radius and luminosity both for Galactic 
globular clusters and for clusters belonging to some other nearby 
Local Group dwarf galaxies. Finally, we discuss what these param- 
eters might be able to tell us about the origins and evolution of 
different types of clusters. 



2 CLUSTER DATA AND CLASSIFICATIONS 

In Table 1 we list all 150 known Galactic globular clusters, includ- 
ing the handful of clusters that are thought to be associated with the 
Sagittarius dwarf galaxy. The structures of most of these objects 
have been measured. For each such cluster we list the half-light ra- 
dius (Rh) and tidal radius (Rt). We have adopted the majority of 
these values from the online database of Harris (1996) (2003 up- 
date). Since this compilation lists Rh and R t as angular diameters, 
we have converted them to parsecs using the cluster distances that 
are also listed by Harris. For a few clusters the database does not 
have entries for Rh and R t . We have been able to fill in these miss- 
ing values for four clusters: IC 1257 (Cote 1999), BH 176, Terzan 
4, and 1636-283 (Webbink 1985). In addition we have, for HP 1, 
adopted the measurements of Ortolani, Bica & Barbuy (1997), who 
used a full colour-magnitude diagram to derive a distance less than 
half that listed by Harris, in good agreement with the measurement 
by Davidge (2000). Also given in Table 1 are the Galactocentric 
distance R gc , the integrated luminosity My, and (where available) 
the ellipticity (e = 1 — b/a, where a is the cluster semi-major 
axis and b the semi-minor axis), again all taken from the Harris 
database. 

Using the half-light radius and integrated luminosity for each 
cluster, we have calculated Ih, the half-light intensity. This is de- 



fined as Ih = „ L „-> , and is measured in units of solar luminosi- 

ZTYKf 
h 

ties per square parsec (Lq pc~ 2 ). It is essentially a measure of the 
mean surface brightness of a cluster within Rh, but in terms of 
absolute units rather than observational units. Because of the evap- 
oration of stars from dynamically evolving clusters, Ih is more sen- 
sitive to evolutionary effects than is Rh. Nevertheless, this param- 
eter turns out to be useful when examining cluster classifications. 
Finally, Table 1 lists the metallicity [Fe/H] and HB morphology in- 
dex of each cluster, along with the formal type classification from 
Mackey & Gilmore (2004). As noted in Section 1, this classifica- 
tion scheme follows closely those of Zinn (1993) and of van den 
Bergh (1993) (see also the references therein). To briefly review 
this scheme, the clusters are grouped by metallicity and HB-type*. 
The metal rich globular clusters with [Fe/H] > —0.8 all have red 
horizontal branches and are spatially confined to the bulge and in- 
ner disk of the Galaxy. Clusters in this sub-group are designated 
bulge/disk (BD) objects. In contrast, the metal-poor clusters (with 
[Fe/H] < —0.8) are seen to have a wide range in HB-type and are 
generally situated in the Galactic halo. 

Most of the inner halo clusters exhibit quite a tight relation- 
ship between HB-type and [Fe/H]. On the other hand many of the 
outer halo clusters have much redder HB-types at given metallicity 
(this is the classical "second-parameter phenomenon", where a pa- 
rameter in addition to metallicity affects HB morphology). Follow- 
ing Zinn (1993), the "old" halo (OH) and "young" halo (YH) sub- 
systems are defined by measuring the offset in HB-type at given 
[Fe/H] of each cluster from the fiducial line which traces the rela- 
tionship between HB-type and [Fe/H] for the inner halo clusters^. 
Objects for which the offset in HB-type is greater than —0.3 (i.e., 
objects which do not lie far to the red of the fiducial line) constitute 
the OH class, while those with offset less than —0.3 (i.e., clusters 
which have red HB-type at given [Fe/H]) form the YH class. 

The majority of the metallicity and HB-type parameters we 
list in Table 1 are taken from the Harris database, except where 
they were updated by Mackey & Gilmore (2004) (one metallic- 
ity value and 28 HB indices). There are five classification types 
listed in Table 1: BD for the bulge/disk clusters (37 objects); OH 
for the old halo clusters (70 objects); YH for the young halo clus- 
ters (30 objects); SG for clusters belonging to the Sagittarius dwarf 
(6 objects); and UN for clusters whose membership type remains 
unknown due to insufficient or unclassifiable data (7 objects). 

In the following two sections, we consider some of the prop- 
erties of these cluster sub-systems which were not considered by 
Mackey & Gilmore, but which nevertheless provide useful insights 
into the origin and evolution of objects in each of the different 
classes of Galactic globular cluster. 



3 PROPERTIES OF THE DIFFERENT GALACTIC 
CLUSTER SUB-SYSTEMS 

It is useful to briefly summarize the results of Mackey & Gilmore 
(2004), who examined some key properties of the Galactic globular 

* HB-type is parametrized by the dimensionless morphology index 

(B+V+fl) (Lee 1990; Lee et aL 1994 )' Where B is the number of HB stars 
which lie to the blue of the instability strip; R is the number of HB stars 
falling to the red of the instability strip; and V is the number of variable HB 
stars 

t Mackey & Gilmore (2004) used the oldest isochrone from the theoretical 
HB models of Rey et al. (2001) as their fiducial 
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Table 1. Data and classifications for Galactic globular clusters. 



Cluster 


-Rgc 


M v 


R h 


Rt 


log h 


e 


[Fe/H] 


HB -index 


Class 


Name 


(kpc) 




(pc) 


(pc) 


(L© pc" 2 ) 


(1 - b/a) 








NGC 104 


7.4 


—9.42 


3.65 


56.10 


3.77 


0.09 


—0.76 


—0.99 


BD 


NGC 288 


12.0 


—6.74 


5.68 


33.12 


2.32 




— 1.24 


0.98 


OH 


NGC 362 


9.4 


—8.41 


2.00 


39.83 


3.89 


0.01 


— 1.16 


—0.87 


YH 


Ndr 12(S1 


18.2 


—7.81 


3.58 


34.73 


3.15 


0.07 


— 1.35 


—0.71 


YH 


Pal. 1 


17.0 


—2.47 


2.16 


28.41 


1.45 


0.22 


—0.60 


— 1.00 


BD 


AM-1 


123.2 


—4.71 


17.73 


68.08 


0.52 




— 1.80 


—0.93 


YH 


LUCU.1 Ui) 


95.2 


—5.14 


10.50 


83.17 


1.14 




— 1.46 


— 1.00 


YH 


Pal. 2 


35.4 


—8.01 


5.38 


54.27 


2.87 


0.05 


— 1.30 


—0.10 


YH 


Mr.r 1851 


16.7 


—8.33 


1.83 


41.18 


3.94 


0.05 


— 1.22 


—0.32 


OH 


Mr.r iqr)4 


18.8 


—7.86 


3.00 


31.30 


3.32 


0.01 


— 1.57 


0.89 


OH 


NGC 2298 


15.7 


—6.30 


2.43 


20.17 


2.88 


0.08 


— 1.85 


0.93 


OH 


NGC 2419 


91.5 


—9.58 


17.88 


214.07 


2.46 


0.03 


—2.12 


0.86 


OH 




41.7 


—5.75 


15.59 


71.37 


1.04 




— 1.30 


— 1.00 


YH 


NGC 2808 


11.1 


—9.39 


2.12 


43.42 


4.23 


0.12 


— 1.15 


—0.49 


OH 


E3 


7.6 


—2.77 


2.58 


13.16 


1.41 




—0.80 




UN 


Pal. 3 


95.9 


—5.70 


17.80 


129.70 


0.91 




— 1.66 


—0.50 


YH 




8.9 


—7.46 


3.90 


41.38 


2.93 


0.12 


— 1.58 


0.08 


YH 


Pal 4 


111.8 


—6.02 


17.15 


105.78 


1.07 




— 1.48 


— 1.00 


YH 


Nnr 4147 

1>IVJ\^ tit/ 


21.3 


—6.16 


2.41 


35.43 


2.83 


0.08 


— 1.83 


0.66 


SG 


Nnr 4372 


7.1 


—7.77 


6.58 


58.75 


2.60 


0.15 


—2.09 


1.00 


OH 


Run 1ftn~ 


18.5 


—6.35 


6.78 


31.02 


2.01 




— 1.67 


—0.82 


YH 


Nnr 4590 

. V I V y\3 


10.1 


—7.35 


4.60 


90.02 


2.74 


0.05 


—2.06 


0.17 


YH 


INVJV^ iOJJ 


7.0 


—8.16 


4.56 


33.75 


3.08 


0.07 


— 1.80 


0.93 


OH 


NGC 5024 


18.3 


—8.70 


5.75 


112.62 


3.09 


0.01 


— 1.99 


0.81 


OH 




16.9 


—6.72 


16.70 


65.21 


1.37 


0.21 


—2.29 


0.52 


YH 


. N V J V J1J7 


6.4 


— 10.29 


6.44 


87.92 


3.63 


0.17 


— 1.62 




UN 


. < V J v JZ / z. 


12.2 


—8.93 


3.39 


115.53 


3.64 


0.04 


— 1.57 


0.08 


YH 


. N v J v JZOU 


8.4 


—8.61 


2.21 


26.75 


3.89 


0.12 


— 1.67 


0.80 


OH 


AM-4 


25.5 


— 1.60 


3.65 


11.48 


0.65 




—2.00 




UN 




16.2 


—6.96 


10.41 


158.36 


1.88 


0.11 


—2.22 


0.58 


YH 


NGC 5634 


21.2 


—7.69 


3.96 


61.28 


3.01 


0.02 


— 1.88 


0.91 


OH 


NGC 5694 


29.1 


—7.81 


3.33 


43.30 


3.21 


0.04 


— 1.86 


1.00 


OH 


IC 4499 


15.7 


—7.33 


8.25 


67.90 


2.23 


0.08 


— 1.60 


0.11 


YH 


NGC 5824 


25.8 


—8.84 


3.35 


144.28 


3.62 


0.03 


— 1.85 


0.79 


OH 


Pal. 5 


18.6 


—5.17 


19.98 


109.87 


0.60 




— 1.41 


—0.40 


YH 


Nnr 5897 


7.3 


—7.21 


7.61 


43.54 


2.25 


0.08 


— 1.80 


0.86 


OH 


Nnr 5904 


6.2 


—8.81 


4.60 


61.96 


3.33 


0.14 


— 1.27 


0.31 


OH 


Nnr 5927 


4.5 


—7.80 


2.54 


36.88 


3.44 


0.04 


—0.37 


— 1.00 


BD 


NCiC S946 


5.8 


—7.20 


2.13 


74.09 


3.35 


0.16 


— 1.38 


0.69 


OH 


BH 176 


9.7 


—4.35 


3.92 


28.16 


1.68 




0.10 


— 1.00 


BD 


1"I\JL. J70U 


4.8 


—8.44 


3.18 


31.83 


3.50 


0.06 


— 1.58 


0.97 


OH 


T vnon 7 


4.2 












—0.62 


— 1.00 


BD 


Pal 14 


69.0 


—4.73 


24.72 


113.07 


0.24 




— 1.52 


— 1.00 


YH 


NGC 6091 


3.8 


—8.23 


1.89 


38.63 


3.87 




— 1.75 


0.93 


OH 


NdC 6191 

, *iv J V U 1 Z. 1 


5.9 


—7.20 


2.34 


20.79 


3.27 




— 1.20 


—0.06 


OH 


Nnr ai 01 

. *i v j v utui 


11.1 


—6.91 


7.61 


32.36 


2.13 


0.05 


— 1.82 


0.84 


OH 


Nnr 6144 


2.6 


—6.75 


4.01 


82.21 


2.62 


0.25 


— 1.75 


1.00 


OH 


Nnr 6139 

lIVJv, yJL^}y 


3.6 


—8.36 


2.41 


25.03 


3.71 


0.05 


— 1.68 


0.91 


OH 


Tpr7an 1 


2.4 


—4.61 


2.84 


12.89 


2.07 




—0.73 


— 1.00 


BD 


Nnr 6171 


3.3 


—7.13 


5.03 


32.47 


2.58 


0.02 


— 1.04 


—0.73 


OH 


1 636-283 

1 UJU Z.O -_J 


2.0 


—3.97 


1.30 


7.84 


2.49 




— 1.50 


—0.40 


YH 


Nnr 6205 


8.7 


—8.70 


3.34 


56.40 


3.56 


0.11 


— 1.54 


0.97 


OH 


NGC 6229 


29.7 


-8.05 


3.27 


47.58 


3.32 


0.05 


-1.43 


0.24 


YH 


NGC 6218 


4.5 


-7.32 


3.08 


25.09 


3.08 


0.04 


-1.48 


0.97 


OH 


NGC 6235 


4.1 


-6.44 


2.79 


25.24 


2.81 


0.13 


-1.40 


0.89 


OH 


NGC 6254 


4.6 


-7.48 


2.32 


27.49 


3.39 




-1.52 


0.98 


OH 


NGC 6256 


1.8 


-6.52 


2.08 


18.55 


3.10 




-0.70 


-1.00 


BD 


Pal. 15 


37.9 


-5.49 


15.70 


64.74 


0.93 




-1.90 


1.00 


OH 


NGC 6266 


1.7 


-9.19 


2.47 


18.00 


4.02 


0.01 


-1.29 


0.32 


OH 


NGC 6273 


1.6 


-9.18 


3.13 


36.27 


3.81 


0.27 


-1.68 


0.96 


OH 


NGC 6284 


7.6 


-7.97 


3.47 


102.72 


3.24 


0.03 


-1.32 


0.88 


OH 


NGC 6287 


2.1 


-7.36 


2.03 


28.43 


3.46 


0.13 


-2.05 


0.98 


OH 
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Table 1. continued... 



Cluster 
Name 


(kpc) 


M 


V 


Rh 

(pc) 


Rt 
(pc) 


log I h 
(Lq pc- 2 ) 


e 

(1 - b/a) 


[Fe/H] 


HB-index 


Class 


NGC 6293 


1.4 


-7 


77 


2 


33 


36.43 


3 


50 


0.03 


-1.92 


0.90 


OH 


NGC 6304 


2.2 


-7 


32 


2 


46 


23.13 


3 


28 


0.02 


-0.59 


-1.00 


BD 


NGC 6316 


3.2 


-8 


35 


2 


27 


18.97 


3 


76 


0.04 


-0.55 


-1.00 


BD 


NGC 6341 


9.6 


-8 


20 


2 


60 


36.18 


3 


58 


0.10 


-2.28 


0.91 


OH 


NGC 6325 


1.1 


-6 


95 


2 


19 


22.08 


3 


23 


0.12 


-1.17 


0.84 


OH 


NGC 6333 


1.7 


-7 


94 


2 


18 


18.75 


3 


63 


0.04 


-1.75 


0.87 


OH 


NGC 6342 


1.7 


-6 


44 


2 


20 


37.17 


3 


02 


0.18 


-0.65 


-1.00 


BD 


NGC 6356 


7.6 


-8 


52 


3 


27 


35.24 


3 


51 


0.03 


-0.50 


-1.00 


BD 


NGC 6355 


1.8 


-8 


08 


2 


40 


41.95 


3 


60 




-1.50 


0.62 


OH 


NGC 6352 


3.3 


-6 


48 


3 


32 


17.43 


2 


68 


0.07 


-0.70 


-1.00 


BD 


IC 1257 


17.9 


-6 


15 


3 


17 


62.57 


2 


59 




-1.70 


1.00 


OH 


Terzan 2 


0.9 


-5 


27 


3 


85 


25.61 


2 


07 




-0.40 


-1.00 


BD 


NGC 6366 


5.0 


-5 


77 


2 


75 


15.92 


2 


56 


0.16 


-0.82 


-0.97 


OH 


Terzan 4 


1.3 


-6 


09 





56 


4.31 


4 


07 




-1.60 


1.00 


OH 


HP 1 


6.1 


-6 


41 


6 


20 


33.71 


2 


12 




-1.55 


0.75 


OH 


NGC 6362 


5.1 


-6 


94 


4 


82 


36.85 


2 


54 


0.07 


-0.95 


-0.58 


OH 


Liller 1 


1.8 


-7 


63 


1 


26 


35.10 


3 


98 




0.22 


-1.00 


BD 


NGC 6380 


3.2 


-7 


46 


2 


33 


37.54 


3 


38 




-0.50 


-1.00 


BD 


Terzan 1 


2.5 


-4 


90 


6 


22 


18.03 


1 


50 




-1.30 


-1.00 


YH 


Ton. 2 


1.4 


-6 


14 


2 


54 


25.33 


2 


78 




-0.50 


-1.00 


BD 


NGC 6388 


3.2 


-9 


42 


1 


95 


18.06 


4 


32 


0.01 


-0.60 


-1.00 


BD 


NGC 6402 


4.1 


-9 


12 


3 


49 


89.92 


3 


69 


0.11 


-1.39 


0.65 


OH 


NGC 6401 


2.7 


-7 


90 


5 


83 


36.96 


2 


76 


0.15 


-0.98 


0.35 


OH 


NGC 6397 


6.0 


-6 


63 


1 


56 


10.58 


3 


40 


0.07 


-1.95 


0.98 


OH 


Pal. 6 


2.2 


-6 


81 


1 


82 


14.35 


3 


33 




-1.09 


-1.00 


YH 


NGC 6426 


14.6 


-6 


69 


5 


78 


79.66 


2 


28 


0.15 


-2.26 


0.58 


YH 


Djorg. 1 


4.1 


-6 


26 


4 


40 


34.98 


2 


35 




-2.00 




UN 


Terzan 5 


2.4 


-7 


87 


2 


49 


39.76 


3 


49 




0.00 


-1.00 


BD 


NGC 6440 


1.3 


-8 


75 


1 


42 


15.42 


4 


33 


0.01 


-0.34 


-1.00 


BD 


NGC 6441 


3.9 


-9 


64 


2 


18 


27.23 


4 


31 


0.02 


-0.53 


-1.00 


BD 


Terzan 6 


1.6 


-7 


67 


1 


22 


48.06 


4 


03 




-0.50 


-1.00 


BD 


NGC 6453 


1.8 


-6 


88 


1 


03 


60.04 


3 


86 


0.09 


-1.53 


0.84 


OH 


UKS 1 


0.8 


-6 


88 


2 


08 


45.29 


3 


25 




-0.50 


-1.00 


BD 


NGC 6496 


4.3 


-7 


23 


6 


26 


17.63 


2 


43 


0.16 


-0.64 


-1.00 


BD 


Terzan 9 


1.6 


-3 


85 


1 


47 


15.54 


2 


34 




-2.00 


0.25 


YH 


Djorg. 2 


1.4 


-6 


98 


1 


62 


20.52 


3 


50 




-0.50 


-1.00 


BD 


NGC 6517 


4.3 


-8 


28 


1 


95 


12.88 


3 


86 


0.06 


-1.37 


0.62 


OH 


Terzan 10 


2.4 


-6 


31 














-0.70 


-1.00 


BD 


NGC 6522 


0.6 


-7 


67 


2 


36 


37.30 


3 


45 


0.06 


-1.44 


0.71 


OH 


NGC 6535 


3.9 


-4 


75 


1 


52 


16.54 


2 


67 


0.08 


-1.80 


1.00 


OH 


NGC 6528 


0.6 


-6 


56 





99 


38.08 


3 


76 


0.11 


-0.04 


-1.00 


BD 


NGC 6539 


3.1 


-8 


30 


4 


08 


52.44 


3 


23 


0.08 


-0.66 


-1.00 


BD 


NGC 6540 


4.4 


-5 


38 





26 


10.21 


4 


45 




-1.20 


0.30 


OH 


NGC 6544 


5.3 


-6 


66 


1 


39 


1.61 


3 


51 


0.22 


-1.56 


1.00 


OH 


NGC 6541 


2.2 


-8 


37 


2 


42 


60.27 


3 


71 


0.12 


-1.83 


1.00 


OH 


2MASS-GC01 


4.5 


















-1.20 




UN 


ESO280-SC06 


14.3 


















-2.00 




UN 


NGC 6553 


2.2 


-7 


77 


2 


71 


14.24 


3 


37 


0.17 


-0.21 


-1.00 


BD 


2MASS-GC02 


4.1 






















UN 


NGC 6558 


1.0 


-6 


46 


3 


47 


22.47 


2 


63 




-1.44 


0.70 


OH 


IC 1276 


3.7 


-6 


67 


3 


69 


32.96 


2 


66 




-0.73 


-1.00 


BD 


Terzan 12 


3.4 


-4 


14 


1 


17 


4.33 


2 


65 




-0.50 


-1.00 


BD 


NGC 6569 


2.9 


-8 


30 


4 


14 


21.63 


3 


22 




-0.86 


—0.82 


OH 


NGC 6584 


7.0 


-7 


68 


3 


12 


36.52 


3 


21 


0.03 


-1.49 


-0.15 


YH 


NGC 6624 


1.2 


-7 


49 


1 


88 


47.22 


3 


58 


0.06 


-0.44 


-1.00 


BD 


NGC 6626 


2.7 


-8 


18 


2 


54 


18.36 


3 


59 


0.16 


-1.45 


0.90 


OH 


NGC 6638 


2.3 


-7 


13 


1 


84 


18.51 


3 


45 


0.01 


-0.99 


-0.30 


OH 


NGC 6637 


1.9 


-7 


64 


2 


20 


22.10 


3 


50 


0.01 


-0.70 


-1.00 


BD 


NGC 6642 


1.7 


-6 


77 


1 


78 


24.61 


3 


34 


0.03 


-1.35 


-0.04 


YH 


NGC 6652 


2.8 


-6 


68 


1 


91 


13.16 


3 


24 


0.20 


-0.96 


-1.00 


OH 


NGC 6656 


4.9 


-8 


50 


3 


03 


26.97 


3 


57 


0.14 


-1.64 


0.91 


OH 


Pal. 8 


5.6 


-5 


52 


2 


14 


50.85 


2 


68 




-0.48 


-1.00 


BD 
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Table 1. continued... 



Cluster 


-Rgc 


M v 


Rh 


Rt 


log I h e [Fe/H] HB-index Class 


Name 


(kpc) 




(PC) 


(PC) 


(L© PC" 2 ) (1 - b/a) 



NGC 6681 


2.1 


-7.11 


2.43 


20.71 


3.20 


0.01 


-1.51 


0.96 


OH 


NGC6712 


3.5 


-7.50 


2.75 


14.93 


3.25 


0.11 


-1.01 


-0.62 


OH 


NGC 6715 


19.2 


-10.01 


3.82 


58.23 


3.97 


0.06 


-1.79 


0.54 


SG 


NGC 6717 


2.4 


-5.66 


1.40 


20.38 


3.10 


0.01 


-1.29 


0.98 


OH 


NGC 6723 


2.6 


-7.84 


4.07 


26.60 


3.05 




-1.12 


-0.08 


OH 


NGC 6749 


5.0 


-6.70 


2.53 


11.97 


3.00 




-1.60 


1.00 


OH 


NGC 6752 


5.2 


-7.73 


2.72 


64.39 


3.35 


0.04 


-1.56 


1.00 


OH 


NGC 6760 


4.8 


-7.86 


4.69 


27.90 


2.93 


0.04 


-0.52 


-1.00 


BD 


NGC 6779 


9.7 


-7.38 


3.41 


25.15 


3.02 


0.03 


-1.94 


0.98 


OH 


Terzan 7 


16.0 


-5.05 


6.55 


49.06 


1.52 




-0.82 


-1.00 


SG 


Pal. 10 


6.4 


-5.79 


1.70 


5.29 


2.98 




-0.10 


-1.00 


BD 


Arp2 


21.4 


-5.29 


15.89 


105.24 


0.84 




-1.84 


0.53 


SG 


NGC 6809 


3.9 


-7.55 


4.46 


25.10 


2.85 


0.02 


-1.81 


0.87 


OH 


Terzan 8 


19.1 


-5.05 


7.56 


30.25 


1.39 




-1.99 


1.00 


SG 


Pal. 11 


7.9 


-6.86 


5.63 


37.13 


2.37 




-0.39 


-1.00 


BD 


NGC 6838 


6.7 


-5.60 


1.92 


10.43 


2.80 




-0.73 


-1.00 


BD 


NGC 6864 


14.6 


-8.55 


2.83 


43.84 


3.65 


0.07 


-1.16 


-0.07 


OH 


NGC 6934 


12.8 


-7.46 


2.74 


38.23 


3.24 


0.01 


-1.54 


0.25 


YH 


NGC 6981 


12.9 


-7.04 


4.35 


45.25 


2.67 


0.02 


-1.40 


0.14 


YH 


NGC 7006 


38.8 


-7.68 


4.59 


76.54 


2.88 


0.01 


-1.63 


-0.28 


YH 


NGC 7078 


10.4 


-9.17 


3.18 


64.42 


3.79 


0.05 


-2.26 


0.67 


YH 


NGC 7089 


10.4 


-9.02 


3.11 


71.75 


3.75 


0.11 


-1.62 


0.92 


OH 


NGC 7099 


7.1 


-7.43 


2.68 


42.68 


3.25 


0.01 


-2.12 


0.89 


OH 


Pal. 12 


15.9 


-4.48 


7.11 


96.78 


1.22 




-0.94 


-1.00 


SG 


Pal. 13 


26.7 


-3.74 


3.45 


23.34 


1.55 




-1.74 


-0.20 


YH 


NGC 7492 


24.9 


-5.77 


9.16 


62.67 


1.51 


0.24 


-1.51 


0.81 


OH 



cluster sub-systems. They found that the majority of globular clus- 
ters associated with nearby dwarf galaxies are, in terms of HB mor- 
phology and core-radius distribution, essentially indistinguishable 
from those in the Galactic YH class. Furthermore, the YH objects 
are observed to be characterized by large, energetic orbits around 
the Galaxy. These orbits are generally of high eccentricity and in- 
termediate inclination, and cover a very large range in orbital angu- 
lar momentum, including highly retrograde orbits. A small fraction 
of OH clusters have similar orbits, but most have much smaller en- 
ergies and eccentricities. Many BD clusters (unsurprisingly) have 
disk-like motions, and/or orbits which are confined to the Galactic 
bulge (e.g., Dinescu et al. (2003)). 

Mackey & Gilmore also examined the distribution of cluster 
ages in each of the sub-systems, using data from the relative age 
study of Salaris & Weiss (2002). These authors determined relative 
ages for 55 globular clusters using several different CMD-based 
techniques. A significant fraction of the YH clusters are ~ 3 Gyr 
younger than the majority of the OH clusters (see Fig. 9 in Mackey 
& Gilmore (2004)). However, several OH clusters are also seen to 
be somewhat younger, while a number of YH clusters appear to be 
just as old as the oldest measured OH clusters. The latter result is 
interesting because it hints explicitly at the possible existence of a 
third parameter, in addition to [Fe/H] and age, which governs HB 
morphology in a cluster. 

Finally, Mackey & Gilmore observed that the YH clusters 
follow a considerably steeper age-metallicity relationship than do 
most of the OH and BD clusters (Mackey & Gilmore, Fig. 10). This 
shows that YH-type clusters apparently underwent much slower 
evolution to higher metallicities than did the OH and BD objects. 

Taken together, the results of Mackey & Gilmore (2004) fully 



support the conclusion of Zinn (1993) that the YH clusters are not 
native members of the Galaxy, but rather are of external origin. Fur- 
thermore, Mackey & Gilmore (2004) showed that some 15 per cent 
of the OH clusters possess core radii, orbital motions, ages and HB 
morphologies which suggest that they too are of external origin. 
The remainder of the OH group, together with the BD system, ap- 
pear to have been formed in the (proto) Galaxy. In the following 
sub-sections, we consider some additional properties of the three 
cluster sub-systems, focusing on structural parameters such as R h 
and ellipticity, as well as on cluster luminosities and surface bright- 
nesses. 



3.1 Distribution with Galactocentric radius 

It has long been known that YH-type clusters typically reside at 
much larger Galactocentric distances than do most OH and BD 
clusters (see e.g., Zinn (1993); van den Bergh (1993); Mackey 
& Gilmore (2004)). This is clearly seen in Fig. 1, which shows 
the cumulative distribution in the number of globulars as a func- 
tion of R gc . Each of the three globular cluster sub-systems are 
quite distinct in terms of their distribution of Galactocentric radii. 
We can quantify the significance of this statement using simple 
Kolmogorov-Smirnov (K-S) tests on each pair of distributions. 
These show that there is only a 2 per cent chance that the BD 
and OH samples were drawn from the same parent population. The 
probability that the YH clusters share the same parent distribution 
as either the OH or BD clusters is even smaller: 0.015 per cent, and 
0.01 per cent, respectively. 

The metal-rich BD clusters are, as expected, the most centrally 
concentrated system, exhibiting a sharp cut-off at R gc = 7.9 kpc. 
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Bulge/Disk Clusters 
Old Halo Clusters 
Young Halo Clusters 
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Galactocentric Radius (R ) (kpc) 

Figure 1. Cumulative distributions in Galactocentric radii R gc for the three 
Galactic globular cluster sub-systems. 




-6 -4 -2 -10 -8 -6 -4 -2 

Integrated Luminosity (M v ) 

Figure 2. Luminosity distributions for the Galactic globular cluster sub- 
systems. 



Only two BD clusters lie outside this radius - BH 176 (i? gc = 9.7 
kpc) and Pal. 1 (R gc = 17 kpc). It is interesting to note that both of 
these clusters have been tentatively associated with the candidate 
dwarf galaxy in Canis Major, and/or the Monoceros tidal stream 
(see e.g., Frinchaboy et al. (2004)). In addition, both van den Bergh 
& Mackey (2004; hereafter Paper I) and Mackey & Gilmore (2004) 
have previously noted that Pal. 1 has an unusually high metallicity 
for its Galactocentric radius. In this respect it strongly resembles 
the Sagittarius clusters Ter. 7 and Pal. 12 (see e.g., Fig. 8 in Mackey 
& Gilmore (2004)). Furthermore, examining Table 1, Pal. 1 and BH 
176 are, respectively, the lowest luminosity (Mv = —2.5) and the 
third lowest luminosity (Mv = —4.4) clusters in the BD ensemble. 
This adds additional weight to the suspicion that these two objects 
might not be native BD clusters 

The OH clusters are considerably more extended in R gc than 
the BD objects, and dominate the inner Galactic halo. However, 
they are not nearly as extended as the YH clusters, which domi- 
nate the outer halo. Only six of the 70 OH clusters (9 per cent) lie 
beyond 7? gc = 20 kpc, and of these, only one (the very unusual 
cluster NGC 2419 - see Section 4.3) has Rgc > 40 kpc. In con- 
trast, one third of the YH system (10 clusters) have R sc > 20 kpc, 
and six of these objects fall outside R gc = 40 kpc. The most dis- 
tant YH cluster is AM-1 at R gc ~ 120 kpc. As noted by Mackey & 
Gilmore (2004), the new classifications presented in Table 1 show 
that the YH clusters also extend into even central regions of the 
Galaxy. These objects are therefore not exclusively a remote-halo 
population. 

3.2 Distribution of cluster luminosities 

In Fig. 2 we plot histograms of the luminosity distributions for clus- 
ters in each of the three sub-systems. From this plot it is seen that 
both the BD and OH systems have peaks near Mv ~ —7.75. 
Furthermore, these two systems contain very few low luminosity 
clusters. Only 11 per cent of BD clusters and 4 per cent of OH 



clusters are fainter than Mv = —5.5. This contrasts with the sit- 
uation for the YH system, which contains a significant population 
of low luminosity objects. It is seen that 27 per cent of the clus- 
ters in this sub-system are fainter than Mv = —5.5. In addition, 
the distribution of YH clusters peaks at somewhat lower luminos- 
ity (Mv ~ — 7) than do the OH and BD distributions. Again, we 
can evaluate the significance of the differences between the three 
distributions using K-S tests. These show that there is only a 5 per 
cent chance that the OH and YH clusters share the same parent 
luminosity distribution. However, the K-S tests are unable to dis- 
criminate between the OH and BD distributions, and the YH and 
BD distributions - the results do not confirm statistically signif- 
icant differences between these parent populations, but they also 
do not confirm statistically significant similarities. This result is at 
least partly due to the relatively small BD and YH sample sizes. As 
discussed above, the primary difference between the distributions is 
the presence (or not) of a small percentage of low luminosity clus- 
ters. Hence, it would be necessary for the sample sizes to be larger 
to confirm (or deny) this difference at a significant level. This is, of 
course, not possible, since we included all known Galactic globular 
clusters in Table 1 . 

We can, however, clearly state that the luminosity distribution 
of YH clusters is significantly broader than that for the OH clusters. 
Combined with our observations from Section 3.1, this result is 
similar to that obtained by Kavelaars & Hanes (1997), who demon- 
strated that outer halo globular clusters (i.e., those metal-poor clus- 
ters with R gc > 8 kpc) have a broader luminosity function than 
the inner halo globular clusters (i.e., those metal-poor clusters with 
R gc < 8 kpc). The present result is, however, distinct from this, 
as we have not used R gc as a discriminator - indeed, Fig. 1 shows 
that there is a significant population of OH clusters with R gc > 8 
kpc, and vice versa for YH clusters. Our observations imply that LF 
broadness is intrinsic to cluster sub-population rather than simply 
location within the Galaxy. 

As discussed by Kavelaars & Hanes (1997), it seems likely 
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that this observed difference in broadness, as well as the apparent 
dearth of low luminosity clusters in the OH and BD systems, can 
(at least in part) be attributed to the destructive influence of tidal 
forces and bulge and disk shocks in the inner Galaxy. A number 
of authors (see e.g., Fall & Zhang (2001)) have shown that grav- 
itational shocks, in addition to cluster relaxation and evaporation, 
are important in shaping the globular cluster mass function. Since 
the YH clusters typically lie at larger Galactocentric radii than do 
OH and BD clusters, they should be much less strongly affected by 
destructive external forces. 

In addition, if the accretion hypothesis is correct, many of 
the YH clusters have spent some fraction of their life outside the 
Galaxy in a more benign tidal environment. Consider, for exam- 
ple, the clusters associated with the Sagittarius dwarf. While it is 
at present very difficult to constrain the early orbital evolution of 
this galaxy, models by Jiang & Binney (2000) have shown that its 
Galactocentric distance at apocentre ~ 1 1 Gyr ago most likely fell 
in the range 60 < R sc < 250 kpc, with the exact value depending 
on its initial mass (specifically, the amount of dark matter it con- 
tained). Models by Helmi & White (2001) and Ibata et al. (2001), 
which assume comparatively small initial masses, show that Sagit- 
tarius may have followed a similar orbit to that presently observed 
(which has apocentric radius ~ 65 kpc, pericentric radius ~ 15 
kpc, and an orbital period ~ 1 Gyr) for the last ~ 10—12 Gyr. 
In each of these scenarios, the Sagittarius dwarf, and its remaining 
globular clusters, spend a significant fraction of their lives away 
from the inner regions of the Galaxy. The observation that four of 
the Sagittarius clusters are still associated with the main body of 
this galaxy provides further evidence that these clusters have not 
suffered prolonged tidal stresses due to the Galaxy. Tidal forces 
strong enough to disrupt a globular cluster are likely to be so strong 
that they would also pull that cluster out of its parent dwarf galaxy. 

Together the factors descibed above may help to explain why 
the YH system has a greater fraction of low luminosity clusters 
than do the OH and BD systems. If destructive effects are indeed 
responsible, then the YH luminosity distribution may well indicate 
that globular clusters are formed with a broad LF. Alternatively, as 
suggested by Kavelaars & Hanes (1997), the shape of the cluster 
LF may be dependent on the environment in which clusters form. 

Finally, it is interesting to examine the combined luminosity 
distribution of globular clusters in the four nearby dwarf galaxies 
studied by Mackey & Gilmore (2004). This distribution is plotted 
in Fig. 3, and bears a strong resemblance to the luminosity dis- 
tribution of the YH clusters. A significant fraction of the external 
globular clusters (25 per cent) are fainter than My — —5.5, and 
in fact this fraction is almost identical to that observed for the YH 
clusters. A K-S test gives a 92 per cent chance that the YH and 
external clusters have the same parent luminosity distribution, but 
only a 3 per cent chance that the OH and external clusters share the 
same parent luminosity distribution. This observation adds another 
example to the list of characteristics in which the YH sub-system 
closely resembles the group of globular clusters associated with 
nearby Local Group dwarf galaxies (Mackey & Gilmore 2004). It 
therefore offers further evidence in favour of the hypothesis that 
much, or all, of the ensemble of YH clusters is of external origin. 

3.3 Distribution of half-light and tidal radii 

Figs. 4 and 5 show the distributions of half-light radii (Rh) and tidal 
radii (Rt), respectively, for the three Galactic cluster sub-systems. 
Just as Mackey & Gilmore (2004) demonstrated for core radii, it is 
clear from these two Figures that the majority of diffuse clusters - 




Integrated Luminosity (M v ) 

Figure 3. Combined luminosity distribution for the globular clusters in the 
LMC, SMC and Fornax and Sagittarius dwarf galaxies. Note how similar 
this distribution is to that for the YH clusters (Fig. 2). 



that is, clusters with large half-light and tidal radii - fall into the YH 
category. Conversely, the large majority of OH and BD objects are 
seen to be compact. Once again, we apply K-S tests to determine 
the significance of the differences between the three distributions. 
These calculations show that there is a less than 1 per cent chance 
that the OH and YH distributions in R^ were drawn from the same 
parent population, and a less than 0.2 per cent chance that the BD 
and YH distributions in R^ share the same parent population. For 
the BD and OH R^ distributions, this probability is 8 per cent. Ex- 
actly similar results are obtained for the R t distributions. 

Given what is known about the distributions in Galactocen- 
tric radius for the three sub-systems (Section 3.1), this result is not 
unexpected when considering R t . The tidal radii of clusters are de- 
termined by how deeply their orbits dive into the the Galactic po- 
tential in which they reside (see e.g., King (1962)). In fact, it has 
been shown that R t is most strongly correlated with the distance of 
a cluster at the pericentre of its orbit (van den Bergh 1994). Hence, 
since the majority of clusters which reside at large R gc belong to 
the YH ensemble, it is to be expected that this group also contains 
the largest fraction of clusters with large R t . By the same token, the 
BD clusters are the most centrally concentrated in 7? gc , and conse- 
quently all the clusters in this sub-sample have comparatively small 
tidal radii. 

The distributions in R^ are somewhat more intriguing. It has 
long been known that the half-light radii of Galactic globular clus- 
ters correlate with their Galactocentric distances, in the sense that 
the most remote clusters are also the largest (e.g., van den Bergh 
& Morbey (1984)). Hence, with only this fact in mind, one might 
expect the YH system to possess most of the clusters with large 
7?h. This is indeed what is observed. However, the fact that the YH 
clusters are strongly suspected of being accreted in a number of 
minor merger events presents a challenge for this explanation. This 
is so because, for clusters which have been more-or-less randomly 
accreted, there is no reason to expect half-light radii to follow the 
same trend with R gc as that which is observed for clusters which 
are native to the Galaxy. We address this puzzle in more detail in 
Section 4.2. 

Irrespective of the above result, the distributions in are in- 
teresting because, as described in Section 1, a cluster's half-light 
radius is stable over many relaxation times. Hence, the observed 
distribution of half-light radii for the YH clusters should, to a large 
extent, reflect the distribution of radii with which these clusters 
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Half-light Radius (R h ) (pc) 

Figure 4. Distributions of half-light radii (Rh) f° r the Galactic globular 
cluster sub-systems. 




Tidal Radius (R t ) (pc) 

Figure 5. Distributions of tidal radii (-Rt) for the Galactic globular cluster 
sub-systems. 

were formed. If this is the case, then the fact that the OH and BD 
systems contain almost no diffuse clusters with Rh > 7 pc al- 
lows two possibilities: either (i) the distributions in Rh for these 
two ensembles also remain mostly unmodified from their original 
states, and the dearth of diffuse objects is due to environmental in- 
fluence on cluster formation (i.e., no diffuse clusters were formed); 
or (ii) the native Galactic cluster population has been strongly mod- 
ified by destructive tidal forces. From the cluster destruction studies 



of Gnedin & Ostriker (1997) and Dinescu, Girard, & van Altena 
(1999) it is clear that destructive influences, such as evaporation 
due to two-body relaxation as well as disk and bulge shocks, must 
be at least partly responsible. These most strongly affect diffuse 
clusters - that is, objects with large Rh and low My. For exam- 
ple, in the formalism of Dinescu et al. (1999) the destruction rate 
due to both bulge and disk shocks is directly proportional to R^, 
and inversely proportional to cluster mass. Furthermore, clusters 
with short orbital periods and small pericentric radii (i.e., many OH 
clusters and most BD clusters) are most strongly affected. In their 
study, Dinescu et al. estimated destruction rates for two of our BD 
clusters, 25 of our OH clusters, and 10 of our YH clusters. They 
provide a list of the nine clusters with the largest destruction rates 
(these objects have expected survival times of ~ 10 Gyr or less). Of 
these, one is a BD cluster, 7 are OH clusters, and one a YH cluster. 
Clearly the native population is most strongly affected. In addition, 
several of the nine have large Rh (e.g., NGC 288, Pal. 5, and NGC 
6144), while two (NGC 6121 and 6397) have been shown to have 
depleted luminosity functions at low stellar masses (Kanatas et al. 
1995; Piotto et al. 1997). This is expected in a cluster if it has been 
subjected to strong external tidal forces. 

Assuming that some fraction of the native Galactic clusters 
has been destroyed as discussed above, it is possible to place con- 
straints on the initial population. Examining Fig. 4, one sees that 33 
per cent of YH clusters have Rh > 7 pc. Mackey & Gilmore (2004) 
concluded that there are presently ~ 100 native Galactic members 
in the BD/OH set (recall that ~ 15 per cent of the OH clusters are 
possibly of external origin). If the initial distribution of cluster sizes 
for native Galactic globulars resembled that presently observed for 
YH objects, then this population constitutes only 67 per cent of the 
original population, since essentially no OH or BD clusters are ob- 
served to have Rh > 7 pc. Hence the original population of native 
Galactic globular clusters may have numbered ~ 150 objects. This 
result is consistent with both the estimate of Mackey & Gilmore 
(2004) who postulated an original Galactic population numbering 
~ 200 objects, and the result of Gnedin & Ostriker (1997) who 
concluded that up to half of the present population of globular clus- 
ters are likely to be destroyed during the next Hubble time. 

It is important to point out that the above argument explicitly 
depends on the assumption that the OH and BD systems originally 
had Rh distributions similar to that which is presently observed for 
YH clusters. Alternatively we again note the possibility that forma- 
tion conditions did not allow the production of diffuse clusters in 
the inner Galaxy. In that case our estimate of ~ 150 original Galac- 
tic globular clusters is an over-estimate. However, not enough in- 
formation about ancient cluster formation conditions is available at 
present to allow this possibility to be confirmed or rejected. 



3.4 Distribution of surface brightnesses 

In Fig. 6 we plot the distributions of Ih, the half-light intensity (es- 
sentially equivalent to the mean surface brightness within _Rh), for 
the three Galactic globular cluster sub-systems. Since this quantity 
is closely related to both Rh and My, it is not surprising to see that 
the vast majority of low surface brightness clusters belong to the 
YH sub-system. The distribution of Ih for this ensemble is quite 
broad with no prominent maximum. In this respect it differs sig- 
nificantly from the much sharper singly-peaked distributions that 
are observed for the BD and OH clusters. A Kolmogorov-Smirnov 
test shows that there is only a 1 per cent chance that the YH and 
BD samples were drawn from the same parent distribution, and a 
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less than 0.1 per cent probability that the YH and OH samples were 
drawn from the same parent distribution. 

The lack of low surface brightness clusters among the OH 
and BD systems results in a scarcity of clusters with /j, < 100 
LQ pc -2 and 7? gc < 15 kpc (as can be seen by inspecting Table 
1). It is also of interest to note that all the nine known Galactic 
globular clusters with 1^ < 20 LQ pc~ 2 are metal-poor objects 
having [Fe/H] < —1.2. As was already pointed out in Sections 
3.2 and 3.3, these effects are likely due to the tidal destruction of 
fragile low surface brightness clusters that ventured too close to the 
massive bulge of the Galaxy. An additional complication is that ob- 
servational selection effects militate against the discovery of low 
surface brightness clusters against the rich foreground of stars in 
the direction of the Galactic bulge. 

It is interesting to note that the location of the peak surface 
brightness for the OH and BD clusters coincides closely with the 
location of the (marginal) surface brightness peak for the YH clus- 
ters, at log 7h ~ 3.25. As previously noted, the YH population has 
apparently not been as strongly affected by tidal forces as have the 
OH and BD populations. This sub-system may therefore provide 
interesting constraints on the initial distribution of cluster surface 
brightnesses (although, as noted previously, ij, is more sensitive to 
cluster evolution than is Rh, because the evaporation of stars, along 
with stellar evolution, causes variations in My). Thirty-seven per 
cent of YH clusters have I h < 100 LQ pc~ 2 , while for the BD 
and OH clusters this fraction is essentially zero. If we again set the 
number of native Galactic members in the BD and OH groups to 
be ~ 100 then, assuming the distribution of surface brightnesses 
for this ensemble used to resemble that now observed for the YH 
clusters, the present population of native Galactic clusters repre- 
sents 63 per cent of the original population, which must therefore 
have numbered ~ 160 clusters. This is consistent with the esti- 
mate already given in the previous Section. It is interesting to note 
that in a thought experiment where we assume all YH clusters with 
Ih < 100 mag pc~ 2 to have been destroyed, the amplitude of the 



surface brightness peak at log 1^ ~ 3.25 rises to ~ 0.32 - which 
is similar to the amplitudes of the peaks that are observed for the 
OH and BD clusters. 



3.5 Cluster ellipticities 

Next we consider the ellipticity distributions of Galactic globular 
clusters, which are plotted in Fig. 7. This Figure shows some inter- 
esting features. In particular, the YH clusters appear on average to 
be somewhat less flattened than the OH clusters, although the small 
YH sample precludes this difference being confirmed as significant 
by a K-S test. Additional measurements will be required to further 
investigate this possible difference. If the average flattening of YH 
clusters does turn out to be less than that of OH clusters, this would 
be unexpected for several reasons. Frenk & Fall (1982) found the 
oldest LMC globular clusters to be slightly more elliptical than typ- 
ical Galactic globular clusters. A similar result was obtained by 
Goodwin (1997). Since a number of other parameters show a strong 
similarity between the YH clusters and the sample of external clus- 
ters (of which LMC objects make up a significant fraction) it is 
perhaps surprising that the YH clusters should be less elliptical on 
average than the OH clusters. Frenk & Fall (1982) also found that 
older LMC clusters are typically less elliptical than younger LMC 
clusters. If this were true in general, one would have expected the 
YH clusters, many of which are ~ 3 Gyr younger than most of the 
OH clusters, to be more elliptical than these Galactic natives. 

Goodwin (1997) has argued that the strength of the tidal field 
in which a cluster evolves is the predominant factor in shaping its 
ellipticity. This is so because a strong tidal field is quite effective 
in reducing the triaxiality of a cluster and in removing any angular 
momentum it may have, thus reducing its ellipticity. This may ex- 
plain why clusters in the LMC are more elliptical, on average, than 
Galactic clusters. However, this prediction appears to be at odds 
with the present observation that the YH clusters, which generally 
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reside at large _R gc and have supposedly spent at least some frac- 
tion of their lifetimes in low-mass galaxies, may typically be less 
flattened than are the OH and BD clusters, which reside in a strong 
tidal field. One possible explanation for this puzzle is that recent 
bulge and disk shocks have re-introduced velocity anisotropy into 
many of the OH and BD clusters (Goodwin 1997). If so, this might 
explain the tail to high ellipticities for the OH clusters, which is not 
present for the YH clusters. 

Finally, it is interesting to note that all three of the lowest 
surface brightness (Jh < 30 Lq pc~ 2 ) clusters for which ellip- 
ticity values have been published are unusually flattened objects 
with 1 - b/a > 0.2. These clusters are NGC 5053, NGC 7492, 
and Pal. 1. Examination of Digitized Sky Survey (DSS) images 
confirms this unusually high ellipticity for NGC 5053 and 7492 
(the DSS image of Pal. 1 is too faint for useful measurements), 
so this result does not seem to be a measurement artifact due to 
low surface brightness. Clearly, more observational data would be 
valuable. Such information could be obtained in the near future 
because there are still a significant number of low surface bright- 
ness clusters for which ellipticity measurements are not yet avail- 
able. In this connection, it is of interest to note that the very faint 
(M v = -3.0lf,;°) and large (i? h = 23 ± 10 pc) globular cluster- 
like object recently discovered in the outer halo of the Galaxy by 
Willman et al. (2004) is of exceedingly low half-light intensity 
(Jh = 0.40lo;3° L pc~ 2 - i.e., logJ h = -0.4+11), and also 
appears to be significantly flattened in its outer regions. 

If the observed ellipticity of the lowest surface brightness 
globulars were due to rotation, then this would hint at a physical 
relation between rotation and surface brightness. Alternatively, it 
is possible that some globular clusters might be embedded in small 
dark matter halos so that the shapes of the lowest surface brightness 
clusters reflect the axis ratios of their dark matter mini-halos. In 
this connection it is noted that Mashchenko & Sills (2004a; 2004b) 
have recently investigated the formation and evolution of globular 
clusters in different types of dark matter mini-halos. They conclude 
that the presence of obvious tidal tails is the only observational ev- 
idence that can reliably rule out significant amounts of dark mat- 
ter in a globular cluster. Furthermore, they write that "cosmologi- 
cal dark matter halos are known to have noticeably non-spherical 
shapes; a stellar cluster relaxing inside such a halo would have 
a spherical distribution in its denser part where stars dominate 
dark matter, and would exhibit isodensity contours of increasingly 
larger eccentricity in its outskirts where DM becomes the dominant 
mass component." It is currently believed (e.g., Bailin & Steinmetz 
(2004)) that massive cosmological dark halos are triaxial with c/a 
ratios of ~ 0.6. However, the numerical resolution on cosmologi- 
cal iV-body simulations is presently insufficient to reliably predict 
the likely flattenings of subhalos with masses that are only a mil- 
lionth of those of the massive halos themselves. Furthermore, the 
shapes of such dark matter subhalos orbiting the Galaxy may have 
been strongly affected by tidal stripping, which rapidly reduces the 
masses of dark matter subhalos once they have been accreted onto 
larger structures (e.g., Gao et al. (2004)). Even so, the case for a de- 
tailed new systematic study of globular cluster ellipticities is strong, 
especially because of the present availability of sensitive wide field 
cameras at a number of large telescopes. Such a study should cer- 
tainly include low surface brightness clusters, because such objects 
are still very poorly studied. It would be particularly interesting to 
examine how ellipticity changes with projected radius, i.e. do clus- 
ters get more or less flattened with increasing projected radius? 
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Figure 8. Plot of log R h vs. My for Galactic globular clusters with R gc > 
15 kpc (filled circles), and those in the LMC (open squares), SMC (star), 
and Fornax and Sagittarius dwarf galaxies (open triangles, and crosses, re- 
spectively). The best fit line, Eq. 1 is marked by the solid line, while the two 
la error lines are marked by the dotted lines. 

4 THE LOG R H VERSUS M v PLANE 

4.1 The link between cluster size and luminosity 

It is informative to consider the correlation between cluster lumi- 
nosity and half-light radius. Several authors have examined this re- 
lationship in the past, arriving at different conclusions. For exam- 
ple, Ostriker & Gnedin (1997) found R h to be strongly correlated 
with luminosity-f for clusters with 5 < R gc < 60 kpc, such that 
Rh oc (L/L©) -0 ' 6 " 5 with correlation coefficient 0.77. However, 
van den Bergh, Morbey & Pazder (1991) found that half-light ra- 
dius also correlates quite strongly with Galactocentric radius R gz 
for clusters with R gc < 40 kpc: Rh oc R g ^' with correlation coef- 
ficient 0.37l ;Jg. They found that if this dependence is taken into 
account, cluster size appears uncorrelated with luminosity - that is, 
a "normalized" radius R^, which is the size a cluster would have 
at given R gc , does not correlate with luminosity. A similar conclu- 
sion was reached by McLaughlin (2000) who found that for the full 
sample of Galactic globular clusters, a weak dependence of Rh on 
luminosity disappears when the variation of Rh on R gc (McLaugh- 
lin found Rh oc i?g C 4 ) is taken into account. However, McLaughlin 
stresses that he cannot reproduce the result of Ostriker & Gnedin 
(1997) for clusters with 5 < R gc < 60 kpc. 

It seems possible that these different conclusions are due, at 
least in part, to different sample selections. McLaughlin (2000) 
considered the full Galactic globular cluster system, and the sample 
of van den Bergh et al. (1991) contained 98 clusters with R gc < 40 
kpc. In contrast, the ensemble of Ostriker & Gnedin (1997) con- 
sisted of 58 non-core-collapsed clusters with 5 < R gc < 60 kpc. 
In Section 3, we discussed the idea that the globular cluster sample 

t assuming a roughly constant M /L ratio for Galactic globular clusters 
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consisting of objects in the central Galaxy has had its properties 
strongly modified by external (tidal) influences. We found support 
for this notion, in agreement with a number of previous authors. 
With this in mind, in the present context it makes sense to con- 
sider a cluster sample which is as free from tidal modification as 
possible. To this end, we consider a Galactic globular cluster sam- 
ple consisting only of objects with R gc > 15 kpc. The majority 
of clusters in this ensemble are YH members. In Paper I we made 
plots of log i? h versus Mv for both this Galactic sample, as well 
as for a sample consisting of globular clusters in the LMC, SMC, 
and Fornax and Sagittarius dwarf galaxies (Figs. 6 & 7 of Paper 
I). Interestingly, both of these plots appear to show a strong in- 
verse correlation between cluster size and luminosity, particularly 
for clusters in the external dwarf galaxies. We reproduce these data 
here, in Fig. 8. 

Apart from four notable outliers, the majority of clusters in 
Fig. 8 exhibit a clear trend in Rh with Mv ■ Three of the deviant 
clusters - Pal. 1, AM-4, and Pal. 13 - are compact but of very low 
luminosity, and lie below and to the right of the majority of the 
data. All three of these objects have been shown to be extremely 
poorly populated globular clusters - AM-4 does not even possess 
any RGB or HB stars (Inman & Carney 1987). Pal. 1 has been 
studied by Rosenberg et al. (1998), who measured an unusually 
flat mass function, very different from those measured for 21 more 
luminous Galactic globular clusters. These authors concluded that 
Pal. 1 has most likely been modified by tidal forces in combina- 
tion with the evaporation of cluster members. Assuming this to be 
the case, they find Pal. 1 must presently be very close to its final 
destruction. Similarly, Siegel et al. (2001) showed that Pal. 13 has 
been strongly affected by tidal forces and is on the threshold of dis- 
solution, a result supported by the unexpectedly high internal ve- 
locity dispersion measured for this object by Cote et al. (2002). We 
note that the inner halo cluster E3 also inhabits this region of the 
log Rh vs. My plane. This cluster is also exceedingly sparse. The 
photometric study by McClure et al. (1985) found a sharp drop- 
off in the cluster luminosity function several magnitudes below the 
main sequence turn off. This supports the suggestion by van den 
Bergh et al. (1980) that E3 has been strongly modified by evapora- 
tive processes. 

Since Pal. 1, AM-4 and Pal. 13 have apparently been strongly 
affected by external destructive forces, we exclude them from our 
present sample and analysis. All three clusters are apparently very 
close to dissolution, so the lower right hand corner of Fig. 8 may 
constitute the ultimate "graveyard" of compact globular clusters. 
This leaves the anomalous object NGC 2419, which is the fourth 
deviant cluster in Fig. 8, and which is discussed in more detail in 
Section 4.3, below. We note that the inclusion or omission of this 
object from the sample does not strongly affect our results (below), 
although the formal errors are reduced if NGC 2419 is not included. 

Applying a least-squares analysis to the sample, one obtains a 
best-fit line: 

log R h = 0.25M V ± 2.5, (1) 

where the formal errors are ±0.03 on the slope and ±0.18 on the 
intercept. This line is plotted in Fig. 8, along with the la error lines. 
These errors represent the case when NGC 2419 is not included in 
the analysis. The trend defined by Eq. 1 is strong, with correlation 
coefficient 0.83. It is important to note however, that the distribu- 
tion of cluster radii in Fig. 8 is dominated by scatter. 

If one replaces the integrated magnitude in Eq. 1 with Mv = 
-2.51og(L/L ) ± Mvq one finds that R h cx (L/Lq)" ' 625 , 
which is almost exactly the result of Ostriker & Gnedin (1997) if 



M/L is assumed to be approximately constant for Galactic globu- 
lar clusters. The physical significance of Eq. 1 is not clear, however. 
One possibility is that the relationship between Rh and Mv has 
been imposed purely through the dynamical evolution of clusters. 
If this is the case, then the majority of evolution must have been in 
Mv as Rh is stable over cluster lifetimes (see Section 1). Mv fades 
as a cluster evolves, due to both stellar evolution and the evapora- 
tion of cluster members. The rate of fading due to stellar evolution 
is regulated by the cluster IMF (because more massive stars die ear- 
lier), while evaporation is caused by both external influences (e.g., 
tidal stripping) and internal relaxation. Loosely bound clusters are 
susceptible to both processes (see e.g., Gnedin & Ostriker (1997)), 
although we note that our sample has been chosen to minimise pos- 
sible modifications due to tidal shocking. In contrast, compact clus- 
ters are significantly more stable, although very compact objects 
suffer accelerated evaporation due to internal relaxation (see e.g., 
Figs. 6 & 7 in Gnedin & Ostriker (1997)). Given the complexity 
of these combined factors, detailed Af-body modelling will be re- 
quired in order to assess how significant evolutionary effects are on 
the Rh vs. Mv diagram, how they vary with cluster concentration, 
and whether they can explain the observed correlation. 

An alternative possibility is that this correlation represents 
some signature of globular cluster formation. If this is so, the de- 
tails of this signature are not immediately evident, as dynamical 
evolution (e.g., fading in Mv) must again be taken into account 
if the presently observed trend is to be traced back to its original 
form. In considering whether Eq. 1 is at least partly due to clus- 
ter formation conditions, it is important to note that a number of 
previous authors (e.g., van den Bergh et al. (1991); McLaughlin 
(2000)) have claimed that a dependence of the half-light radius R h 
on Galactocentric distance R gc might contribute to the observed 
correlation between cluster size and luminosity for Galactic glob- 
ular clusters. A physical interpretation of this argument is that the 
properties of the Galactic globular cluster system were set by the 
cluster formation process, which is environment-dependent - that 
is, the parameters with which a cluster was formed were modulated 
by forces that depend on Galactocentric radius (see e.g., McLaugh- 
lin (2000) for discussion). The basis for this idea can be seen in 
Figure 9, which shows a plot of Rh as a function of R gc for the 
full sample of Galactic globular clusters. There is a clear trend of 
increasing Rh with increasing R gc . A formal fit to these data yields 

log R h = 0.42 log R gc + 0.19, (2) 

where the errors are ±0.04 on the both the slope and intercept, 
and the correlation coefficient is 0.65. This result is consistent with 
those found previously by both McLaughlin (2000) and van den 
Bergh et al. (1991). 

In our present sample however, the extent to which a correla- 
tion between Rh and Galactocentric radius can drive that between 
Rh and Mv is not clear. All our Galactic clusters are in the outer 
halo, with R gc > 15 kpc. Fig. 9 shows that these objects have a 
very large scatter in R h , and furthermore that diffuse clusters with 
R gc ~ 15 kpc do not have systematically smaller Rh than diffuse 
clusters with i? gc ~ 100 kpc. Any apparent trend in Rh with R gc 
for these outer halo clusters is most likely caused by the surpris- 
ing lack of compact clusters at very large Galactocentric radii (see 
Section 4.2 below). In addition, the majority of the sample with 
Rgc > 15 kpc are YH members, and as such are very possibly ac- 
creted objects. In this scenario, such clusters would not be expected 
to exhibit any trend between R h and R gc amongst themselves. We 
also note that if the data in Fig. 8 are split into the Galactic and 
"external" clusters, then fitting these two groups separately repro- 
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Figure 9. Plot of log _R h vs. log R gc for all Galactic globular clusters. The 
BD clusters are marked by open triangles; the OH clusters by crosses; the 
YH clusters by filled circles; the Sagittarius clusters by open circles; and 
the unclassified (UN) clusters by open squares. The vertical dotted line in- 
dicates R gc = 15 kpc. 

duces Eq. 1 for each to within la, albeit both with lower signif- 
icance (due to the smaller number of points in each group). It is 
curious that the globular clusters in the LMC, SMC, and Fornax 
and Sagittarius dSph galaxies follow such a similar trend in Fig. 8 
to the outer halo Galactic clusters, even though they are members 
of significantly different galaxies and have evolved in different en- 
vironments. Even so, it is observed that in the LMC, the most dif- 
fuse clusters (Reticulum, ESO121-SC03, NGC 1841, 2257) all lie 
at large galactocentric radii, as does the very extended cluster 1 in 
the Fornax dSph, so there seems to be at least a weak dependence 
of Tih on radius in these galaxies. 

Ultimately, the most important conclusion to draw from Fig- 
ures 8 and 9 is that the Galactic globular clusters with large Rh 
are both faint, and occur at large Galactocentric radii. With present 
data, it is not possible to identify for certain the origin of the ob- 
served trend in Rh with My. Although there is some evidence that 
a correlation between R^ and R gc could be at least partly respon- 
sible, at present an intrinsic relationship between Rh and My, per- 
haps a relic of the cluster formation process, cannot and should not 
be ruled out. This puzzle seems well suited for study using large- 
scale realistic iV-body simulations (see e.g., Aarseth (1999)), which 
with the advent of powerful dedicated hardware (e.g., GRAPE-6 - 
see Makino et al. (2003)) are now nearly capable of direct globular 
cluster modelling. 

4.2 Where are the compact outer halo globular clusters? 

As noted above, examining Fig. 9 in more detail reveals a surpris- 
ing result: no compact YH clusters occur at R gc > 40 kpc. This 
observation is intriguing because the presently accepted paradigm 
is that the young halo objects have been accreted by the Galaxy 
via minor mergers. Their Rh values should therefore be indepen- 



Figure 10. Plot of log I h vs. log R gc for all Galactic globular clusters. The 
point styles are as in Fig. 9. 

dent of i? gc . This is so because the chaotic acquisition of clusters 
is expected to result in a wide spread of cluster sizes across a broad 
range of R gc values. Indeed, this is exactly what is seen at interme- 
diate Galactocentric radii (10 < R gc < 40 kpc). The dispersion in 
Rh for the YH clusters at such radii is considerably larger than that 
for the clusters that are intrinsically associated with the Milky Way, 
and in fact closely matches the spread in R h among the Sagittarius 
clusters. At very small R gc we only really see compact YH clusters. 
However, recalling the apparent influence tidal forces have had on 
the OH and BD samples, this is not unexpected because extended 
clusters are especially vulnerable to rapid tidal destruction in the 
inner Galaxy. 

It is therefore the handful of YH clusters at very large Galacto- 
centric radii § (_R gc > 40 kpc) which are the most intriguing. A sim- 
ple examination of Fig. 4 shows that, among a randomly selected 
sample of six YH clusters, one should expect approximately three 
to have Rh < 5 pc, or more than four to have Rh < 10 pc. How- 
ever, Fig. 9 shows that there are no compact clusters among the six 
observed at very large Galactocentric radii. If the very outer young 
halo clusters originally belonged to now-defunct dwarf galaxies, as 
the presently accepted paradigm would seem to suggest, then there 
is no clear reason why their dispersion in Rh should not be similar 
to that observed for YH clusters with 10 < R gc < 40 kpc. 

Inspection of Table 1 reveals an additional interesting fact. All 
of these six outer YH clusters are of relatively low luminosity (in- 
deed, this result might be expected given Fig. 8 and Eq. 1). The sur- 
face brightness Ih is derived from the parameters Rh and My, and 
is therefore useful in the present analysis. Fig. 10 shows a plot of I h 
as a function of R gc for all Galactic globular clusters. In the inner 
Galaxy, there is no significant trend in Ih with R gc , but (as dis- 
cussed in Section 3.4) there is a marked lack of low surface bright- 
ness clusters - which are the most difficult to find against the back- 

§ AM-1, Eridanus, Pyxis, Pal. 3, Pal. 4, and Pal. 14 
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drop of the Galactic bulge, and which are most vulnerable to tidal 
disruption. At intermediate R gc , the YH clusters are seen to have 
a significantly larger spread in Jh than do the native Galactic clus- 
ters. In fact, the dispersion in I h observed at these Galactocentric 
radii for the YH clusters matches closely the spread in Jh exhibited 
by the Sagittarius objects. At 7? gc > 40 kpc only very low surface 
brightness clusters are found, i.e. only low luminosity clusters with 
large Rh- Apart from NGC 2419, which is also interesting for other 
reasons (see Section 4.3, below), there are no globular clusters of 
even average surface brightness in the outer Galactic halo. 

We can think of several possible scenarios that might account 
for this lack of compact/luminous clusters at very large R gc : 

(1) Clusters with R gc > 40 kpc once belonged to one, or more 
likely several, now-destroyed parent dwarf galaxies which disinte- 
grated at large R gc and which only possessed diffuse clusters at the 
time of merger. 

(2) Clusters with R gc > 40 kpc were outlying objects in one, or 
more likely several, parent dwarf galaxies, and were lost very early 
on during merger events which subsequently proceeded at smaller 

R gc . 

(3) Clusters with R gc > 40 kpc never belonged to a parent 
dwarf galaxy and were formed as isolated objects at large Galac- 
tocentric radii. 

Since little is known about the properties of galaxies which 
might have merged with ours at early epochs, it is not possible to 
definitively rule out any of these options. If these galaxies were 
similar to the Fornax and Sagittarius dSph galaxies observed today, 
then options 1 and 2 seem implausible. One could perhaps imag- 
ine a scenario in which any compact cluster members of a dwarf 
galaxy were formed close to the centre of the parent galaxy and 
soon merged into its core, leaving only diffuse clusters. However, 
both the Fornax and Sagittarius dwarf galaxies contain clusters with 
a wide range of Rh, so, unless the hypothetical parent dwarf galax- 
ies of the outer halo globulars were radically different to Fornax 
and Sagittarius, option 1 appears unlikely. In addition, Fornax is at 
a distance of R gc ~ 140 kpc and shows no signs of being strongly 
tidally influenced by the Galaxy, let alone disrupted. It also pos- 
sesses two outlying clusters (Fornax 1 and 5) which have not been 
removed. Furthermore, the Sagittarius dwarf is much closer to the 
Galactic centre (R gc ~ 20 kpc), and still apparently possesses four 
clusters, of which two are diffuse objects (Terzan 8 and Arp 2). 
Hence, option 2 also appears unlikely in this scenario. 

Having outlined this argument, we point out that we cannot 
be certain that the galaxies which merged with ours at early epochs 
were in any way similar to Fornax and Sagittarius. Indeed, studies 
of the chemical abundance patterns in halo stars and in local dSph 
galaxies suggest that perhaps they were not (see Mackey & Gilmore 
(2004) for a brief summary). It is possible that the hypothetical 
parent dwarf galaxies of the outer YH clusters, along with their 
member clusters, were systematically more diffuse than Fornax and 
Sagittarius and their clusters. If this were the case, it might explain 
why these diffuse galaxies are now destroyed whereas Fornax and 
Sagittarius are still (mostly) intact. In this scenario, we cannot rule 
out options 1 and 2. 

Option 3 also seems plausible, as the outermost YH clusters 
might help form the upper end of a trend in increasing Rh with 
increasing R gc if considered together with the BD and OH clus- 
ters (Fig. 9). However, at least four of these outer young halo clus- 
ters are believed to be several Gyr younger than the oldest Galactic 
globulars (e.g., Stetson et al. (1999); Sarajedini (1997)). This raises 
the question: why (and how) did cluster formation occur at such re- 
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Figure 11. Plot of log R h vs. My for all Galactic globular clusters. The 
point styles are again as in Fig. 9. The giant M31 globular cluster Gl is 
marked by two linked stars - these represent the two discrepant measure- 
ments of Rh for this cluster. The upper point is the measurement of Meylan 
et al. (2001), while the lower point is that of Barmby et al. (2002). 

mote locations after it occurred deeper within the Galactic potential 
well. Furthermore, it is not clear why (or how) such cluster forma- 
tion would produce only diffuse, low luminosity clusters, and why 
these objects are similar in so many ways to those found in local 
dwarf galaxies. 

The outermost clusters are clearly very important tracers of 
how remote parts of the Galactic halo were formed or assembled, 
as well as the properties of the systems in which they might have 
formed. It is therefore clearly of great interest to study them in con- 
siderably more detail than has been achieved to date. 

4.3 The structures of very luminous clusters 

Figure 1 1 shows a plot of log Rh vs. My for all Galactic globular 
clusters. The most important thing to note from this Figure is that 
the upper left envelope of the main distribution of globular clusters 
has a rather sharp edge. Only a very few clusters lie a significant 
distance above this upper edge. If we shift the line defined by Eq. 
1 upwards so that it has intercept 2.95 (as plotted on Fig. 1 1) then 
one finds that the line 

log R h = 0.25M V + 2.95 (3) 

closely describes the upper envelope of the cluster distribution. 
Only three Galactic globular clusters lie above this line: u Cen- 
tauri, M54, and NGC 2419. A number of authors, such as Bekki & 
Freeman (2003), Tsuchiya et al. (2003), and Ideta & Makino (2004) 
(and references therein), have demonstrated that the observed prop- 
erties of uj Centauri are consistent with it being the stripped core of 
a former dwarf spheroidal or dwarf elliptical galaxy. Furthermore, 
Ibata, Gilmore & Irwin (1994; 1995) have shown that M54 (= NGC 
6715) is associated with the central region of the Sagittarius dwarf 
galaxy, and some authors (see e.g. Layden & Sarajedini (2000)) 
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have argued that M54 may actually be the core of this galaxy. Fi- 
nally, in Paper I we suggested that NGC 2419 might also be the re- 
maining core of a defunct dwarf spheroidal-type galaxy. In this con- 
nection it would clearly be of great interest to check if NGC 2419, 
like u) Centauri (e.g., Bedin et al. (2004); Hilker et al. (2004); and 
references therein), also shows evidence for age/metallicity sub- 
structure in its color-magnitude diagram. However, evidence for 
such substructure in NGC 2419 might be difficult to obtain because 
the mean metallicity ({[Fe/H]) = -2.12) ofNGC2419is so much 
lower than that of w Centauri (([Fe/H]) = —1.62). As a result the 
expected range in metallicities for the population components of 
NGC 2419 is expected to be quite small. NGC 2419 is also sig- 
nificantly more distant than u> Centauri (R@ — 84.2 kpc), further 
complicating any observations. 

It is intriguing that uj Centauri, M54, and NGC 2419, which 
have all independently been previously suggested as the cores of 
now-disrupted dwarf galaxies, should constitute a group of objects 
which have such strikingly different structures than other lumi- 
nous Galactic clusters. This observation suggests a similar origin 
for these three objects, since each has clearly evolved under quite 
different circumstances and in different environments. Inspection 
of Table 1 shows that three out of the four Galactic globular clus- 
ters brighter than My = — 9.5 fall above Eq. 3. The lone exception 
is NGC 6441 (M v = -9.64, R h = 2.2 pc). The high metallicity 
([Fe/H] = -0.53) of this cluster, which is located at R gc = 3.9 
kpc, militates in favor of its assignment to the population associated 
with the Galactic bulge or inner disk. It is noted in passing that the 
apparently normal globular NGC 5053 also falls very close to (but 
still marginally below) the line defined by Eq. 3. 

Since the majority of the most luminous Galactic globular 
clusters are now suspected of being the cores of former dwarf 
galaxies it seems likely that the same might also be true in the glob- 
ular cluster system associated with M31. This suspicion is strength- 
ened by the observation that the cluster M31 Mayall No. II (= Gl), 
which may be the most luminous globular cluster associated with 
the Andromeda galaxy, has Mv = —10.94 and Rh = 14 pc (Mey- 
lan et al. 2001), although Barmby, Holland, & Huchra (2002) find 
Rh = 3.1 pc (using the distance of Meylan et al. (2001)). Adopt- 
ing either of these Rh measurements, Gl also falls above the line 
defined by Eq. 3, i.e. in the same region as u> Cen, M54, and NGC 
2419. In this connection, we note that Meylan et al. (2001) found 
their HST/WFPC2 photometry of Gl to be consistent with an in- 
ternal metallicity spread of ~ 0.5 dex. Furthermore Bekki & Chiba 
(2004) have recently argued that the observed properties of Gl can 
be explained if this object is the tidally stripped nucleus of a former 
dwarf elliptical (dE,N) galaxy. 

Finally, Harris et al. (2002) constructed a catalogue of struc- 
tural measurements for globular clusters in NGC 5 128 (= Centaurus 
A). The 14 most luminous of these have recently been studied spec- 
troscopically by Martini & Ho (2004), who hypothesize that some 
of these clusters are actually the stripped nuclei of dwarf galaxies, 
based on their large masses and the possible detection of extra-tidal 
light by Harris et al. (2002) for several. On the assumption that our 
line Eq. 3 provides a reasonable means of discriminating objects 
which may be the cores of former dwarf galaxies, the speculation 
of Martini & Ho (2004) is supported by the observation that all 
14 of the clusters in their sample lie above and to the left of this 
line. In other words, all 14 have structures and luminosities more 
similar to those of to Cen, M54, NGC 2419, and Gl than to those 
of the majority of "ordinary" globular clusters associated with our 
Galaxy. Harris, Harris & Geisler (2004) have recently found that 
NGC 5128 contains 980 ± 120 globulars, so that potential dwarf 



galaxy remnants represent at least ~ 1.4 per cent of the clusters 
associated with this galaxy. 



5 CONCLUSIONS 

In terms of their observed characteristics (metallicities, HB mor- 
phologies, structures, shapes, luminosities, and surface bright- 
nesses) the Galactic globular clusters, considered as a whole, con- 
stitute a relatively inhomogeneous class of objects. From a physi- 
cal point of view, it therefore makes sense to try to identify sub- 
groups made up of clusters with similar characteristics. One of 
the most successful classification schemes for the Galactic glob- 
ular clusters is that introduced by Zinn (1993) and van den Bergh 
(1993) (see also the references therein), whereby clusters are segre- 
gated according to their metallicity and HB morphology. Recently 
this scheme has been updated to cover the entire Galactic globular 
cluster system by Mackey & Gilmore (2004). 

In the present paper, we have presented the complete list of 
cluster classifications of Mackey & Gilmore (2004). Several objects 
in the Galactic halo possess properties which are consistent with the 
hypothesis that they are the stripped cores of former dwarf galax- 
ies (see below). The remaining majority of "normal" globular clus- 
ters fall into three major categories - bulge/disk (BD), "old" halo 
(OH), and "young" halo (YH) clusters. As recently re-emphasized 
by Mackey & Gilmore (2004) (see also Section 3, above), the ma- 
jority of BD and OH objects seem to have formed in the (proto) 
Galaxy, while the YH objects possess characteristics which indicate 
that they most likely formed in extra-Galactic systems and were ac- 
creted into the Galactic halo during subsequent merger events. 

We have examined several of the properties of these three sub- 
systems which were not considered by Mackey & Gilmore (2004). 
The BD, OH, and YH clusters typically occupy different regions of 
the Galaxy. The metal-rich BD objects are the most centrally con- 
centrated, exhibiting a sharp cut-off in their Galactocentric radial 
distribution at R gc — 7.9 kpc. In contrast, the OH clusters dom- 
inate the inner halo, with the majority residing within R gc — 20 
kpc, while the YH clusters typically occupy the outer halo, with 
more than half lying at R gc > 15 kpc. 

In terms of cluster structures, luminosities, surface bright- 
nesses, and ellipticities, the YH clusters are clearly distinct from 
the OH and BD objects. The OH and BD populations have appar- 
ently been strongly modified by destructive tidal forces and bulge 
and disk shocks, so that these sub-systems possess few extended 
clusters, few low luminosity clusters, and few low surface bright- 
ness clusters. In comparison, the YH group contains a significant 
fraction of extended clusters with half-light radii Rh > 7 pc, and 
low luminosity clusters with Mv > —5.5. The luminosity distri- 
bution of YH clusters matches closely that observed for clusters in 
four nearby dwarf galaxies (the LMC, SMC, and Fornax and Sagit- 
tarius dwarf spheroidals), and as such, offers further evidence for 
the external origin of much or all of the YH ensemble. By exam- 
ining how the observed distributions of half-light radii and mean 
surface brightnesses for the YH clusters would have to be modified 
to match those observed for native Galactic globulars, we estimate 
that the original population of Galactic globular clusters may have 
numbered ~ 150 — 160 objects. This result is consistent with sev- 
eral previous estimates in the literature. 

The distribution of cluster ellipticities shows that the YH clus- 
ters are, on average, less flattened than are the OH and BD clus- 
ters. A possible explanation of this result is that recent bulge and 
disk shocks have introduced significant velocity anisotropy (and 
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hence ellipticity) into many BD and OH clusters (Goodwin 1997). 
We also note that all three of the lowest surface brightness clusters 
for which ellipticity values have been published are unusually flat- 
tened objects, as is the newly discovered very low surface bright- 
ness object of Willman et al. (2004). This may indicate an intrin- 
sic relationship between cluster rotation and surface brightness. An 
alternative, and perhaps more attractive, explanation invokes the 
influence of dark matter. In the globulars of lowest surface den- 
sity the structure of the potential energy well might be dominated 
by a triaxial dark matter mini-halo (Mashchenko & Sills 2004a; 
Mashchenko & Sills 2004b), rather than by the baryonic material 
associated with Population II stars. The flattening of the isophotes 
of low surface brightness globular clusters might therefore provide 
information on their (still hypothetical) dark matter mini-halos. We 
briefly made the case for a comprehensive study of Galactic globu- 
lar cluster ellipticities. 

We have also investigated the relationship between cluster lu- 
minosity and half-light radius. A strong correlation between the 
two is observed for clusters in the outer Galactic halo and also 
among clusters in the LMC, SMC, and Fornax and Sagittarius 
dwarf galaxies. Several authors have previously asserted that any 
correlation between i?h and Mv for Galactic globular clusters can 
be explained by the fact that also appears to correlate with 
Galactocentric radius. However, this cannot explain why the outer 
halo Galactic clusters, which are mostly YH objects, follow the 
trend in R h with My. If these objects are accreted, no trend in R h 
with R gc is expected. 

An additional puzzling observation is the lack of compact 
globular clusters in the remote halo. We observe six YH clusters 
with R gc > 40 kpc; however all are diffuse and of low luminosity. 
Based on the observed distributions of Rh and Mv for YH clus- 
ters, one might have expected several of the six to be compact and 
of greater than average luminosity. We discussed several scenarios 
which might explain the dearth of compact clusters with very large 
R gc ; however none appear particularly attractive. 

Finally, if all the Galactic globular clusters are plotted on the 
log R h vs. Mv plane, several clusters stand well clear of the main 
body of objects. At very low luminosities we observe a number of 
compact objects (Pal. 1, Pal. 13, E3 and AM-4) which have appar- 
ently been strongly affected by shrinkage due to cooling and evapo- 
ration of cluster stars, and stripping by Galactic tides. The most in- 
teresting result, however, concerns the upper envelope of the main 
distribution of clusters, which is sharply defined. Only three Galac- 
tic globulars lie significantly above a line (Eq. 3) describing this 
upper envelope - to Cen, M54, and NGC 2419. All three of these 
objects have previously, and independently, been hypothesized as 
being the stripped cores of former dwarf spheroidal-type galaxies. 
This suggests that the position of a cluster in the log i?h vs. Mv 
plane provides a powerful tool for the segregation between "nor- 
mal" globular clusters and those which might previously have been 
associated with the cores of dwarf galaxies. This suspicion is rein- 
forced by the additional observation that the massive cluster Gl in 
M3 1 , which has been previously suggested to be the remaining core 
of a nucleated dwarf elliptical galaxy, also lies above Eq. 3, as do 
14 of the most luminous clusters in NGC 5128, which have again 
been postulated as constituting the remnant cores of now defunct 
dwarf galaxies. 
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